ROM2F/2004/13 



Layered Higgs Phase as a Possible Field Localisation on a Brane 

P. Dimopoulos'^")* and K. Farakos^*^^ 

(«) INFN-Rome2 Universita di Roma 'Tor Vergata' 
Dipartimento di Fisica 1-00133, Rome, Italy 
(b) Physics Department, National Technical University 
15780 Zografou Campus, Athens, Greece 

Abstract 

So far it has been found by using lattice techniques that in the anisotropic five- 
dimensional Abelian Higgs model, a layered Higgs phase exists in addition to the 
expected five-dimensional one. The exploration of the phase diagram has shown that 
the two Higgs phases are separated by a phase transition from the confining phase. This 
transition is known to be first order. In this paper we explore the possibility of finding 
a second order transition point in the critical line which separates the first order phase 
transition from the crossover region. This is shown to be the case only for the four- 
dimensional Higgs layered phase whilst the phase transition to the five-dimensional 
broken phase remains first order. The layered phase serves as the possible realisa- 
tion of four-dimensional spacetime dynamics which is embedded in a five-dimensional 
spacetime. These results are due to gauge and scalar field localisation by confining 
interactions along the extra fifth direction. 
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1 Introduction - Motivation 



Since the mid eighties lattice gauge models with anisotropic couplings defined in higher 
D-dimensional spaces have been proposed. These models may exhibit, through a phase 
transition, a phase which is coulombic in (D-1) dimensions and shows confinement along 
the remaining dimension. In fact, this was the result of Fu and Nielsen using mean field 
techniques in a five-dimensional pure U(l) gauge theory with anisotropic couplings pp. This 
new phase was called layered. 

The Monte-Carlo analysis which followed j2] supported the mean field results and helped 
to get a more precise picture of the phase diagram [^. Also in jl] the orders of the phase 
transitions have been analysed 

In addition, as it may have been expected, the consideration of the interaction with 
a scalar particle leads to a richer phase diagram. Actually, the exploration of the phase 
diagram of the model for various sets of lattice parameters values provides strong evidence 
that the layer phase is stable and it appears either in a Higgs phase for the U(l) case [HIH], 
or in a Coulomb phase for a SU(2) adjoint Higgs model § fn\ . 

Since gauge theories defined on a D > 4 spacetime are known to be non-renormalisable 
an explicit cut-off A has to be introduced J2] ■ Therefore the theory is to be considered as an 
effective theory which emerges from a more fundamental renormalisable theory (for example 
the string theory). For the U(l) gauge field the introduction of the cut-off A leads to the 
admission of the strong coupling phase to be the interesting phase for the five-dimensional 
theory. As a consequence the lattice methods have to be used as the unavoidable non- 
perturbative tool for the study of the system. 

Up to now the Monte-Carlo results show that the transition between the five dimensional 
strong coupling phase and the layered Higgs phase is first order. A multilayer structure arises 
which supports the idea of the confinement along the extra dimension [HI E]- A crucial 
question may arise: is there any possibility for this phase transition to be of second order? 
We work on this possibility and we look for a second order ending point along the first order 
critical line ^. This would give evidence for the layer mechanism to be more realistic and 
useful in scenarios concerning the localisation of the fields on the four-dimensional subspace. 

■''It has to be noticed that for non-Abelian gauge theories the layer phase exists in six dimensions 
For the lattice realisation of the four-dimensional confining phase in a five-dimensional non-abelian gauge 
theory in the context of a compactified extra dimension the reader may refer to |S1 E) • 

''Recently a paper appeared [HI which presents a non-perturbative study of the Dvali-Shifman mechanism 
[in] of the gauge localisation on a brane. For that reason a SU(2) gauge theory with an adjoint scalar, whose 
mass parameter is space dependent, is employed in 3D. 

similar behaviour has been seen in U(l)-Higgs model in 4D ^^li in SU(2)-Higgs model in 3D ^31 and 
in SU(2) adjoint Higgs model in 3D |15| . 
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Before proceeding to the lattice model let us present the action of the U(l)-Higgs model 
in five dimensions which in principle could inspire the lattice action used in the sequel for 
the numerical simulation. 

We assume a five dimensional anti de Sitter space (AdSs) with one warped extra di- 
mension. In general the metric reads: 

ds'^ = a^iz)[dxl-dx'^]-dz'^ (1.1) 

We consider 77^,^ to be the four dimensional Minkowski metric and a{z) the warp factor. We 
do not need to define explicitly the form of the warp factor. We only require that it goes to 
zero as 2; — s> 00 (^Hj,^!!, ^Hj, C^])- Hence the five-dimensional metric is written: 

^7^^ = (^^'^^ -1) (1-2) 

We consider now that in such a space we define a five dimensional Abelian Higgs model, 
the action of which reads: 

S = -^gauge + -^scalar 

^ J d'x^ FMNFKLQ^'^g'"' + j d'x^ [DM<^*D^<!>g^''' - V{^)] 



d'^xdz 



- — F F xn'^'^n''^ - ^ F ^F ^n^" 



+ 



^ j^.,. [aH.)D^^-n.^,r - .H.)D.^-D.^ - .W)] ,1.3) 

We note that the upper case indices refer to the 5-D space, M,N,K,L = 0, ...,4 and the 
lower case Greek ones to the 4-D space i.e. fi, u, n, X = 0, 3. It is obvious that the scalar 
field $ depends on the five dimensional space {x, z). Then we use the rescaling: a{z)(^ = (p 
for the scalar field. In the rather general case where the quartic scalar potential is considered, 
the scalar action takes the form: 

•^scalar = / ^^^^z [D^^*D^ip - a^{z)D,^*D,^ - M{zf^*^ - (1.4) 

where M\z) = c?{z)m^ + \a'{z)f + ^[^^(z)]" H. 

II Assuming that rr? < on the brane (z = 0), we note that depending from the exact form of the warp 
factor the mass term may turn to be positive after a certain distance or at least tends to zero asymptotically 
along the transverse direction. So we meet the situation of two degenarate minima near the brane and only 
one minimum far away from it. 
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It is a trivial matter for the action to be analytically continued to the Euclidean space 
from which the lattice action can be defined after following the usual methods for discretiza- 
tion. Therefore we take: 



Sl - Sgauge + 'S'scalar 

= /^^E E (l-cost/^.(x)) + 5^ 5^ (l-cosf/^5(x)) 

X l<fi<i/<4 X 1</J<4 

+ ^ l^Lix) - Uf,{x)LpLix + afi)Y[LpLix) - Uf,{x)LpLix + afi)] 

X 1<M<4 

X 

+ Y.'^^L^L{^)^L{x)+PR{^l{x)VL{x)fl (1.5) 

X 

We denote by (^l^x) the lattice scalar field and 

Uf,y{x) = U^{x)Uy{x + afi)Ul{x + ai))Ul{x) 

U^5{x) = U^{x)U^{x + afi)Ul{x + ab)Ul{x) (1.6) 

are the plaquettes on the four-dimensional space and along the fifth direction respectively 
The U's are the links for the gauge field on the lattice **. They are explicitly given by: 
Um = 6^°"^'^' (with M = 1, 5). The primed couplings refer to the interactions along the 
extra dimension. Moreover as it can be noticed from the corresponding continuous action, 
the couplings obey certain relationships, which depend on the warp factor Hence we have: 

0, = a\xT)^, /3; = a2(xT)/3., A = ^ (1.7) 

a^M\xT) = ^ml (1.8) 

Therefore due to the assumed form for the warp factor the interactions for both the gauge 
and scalar fields are strongly coupled along the extra direction. 

Since a brane is defined as any three dimensional submanifold to which ordinary matter is 
trapped so that it can not escape to the bulk, a possible realisation of the trapping 
mechanism is to assume the existence of confinement along the extra dimension. On the 
lattice this situation can be realised using a lattice model with anisotropic couplings. This 



**Notice also that here we use the symbol x for the whole discretised five-dimensional space. The extra 
direction now is denoted by xt- 

ttPor the transition from the continous to the lattice action we have assumed the following rescaling for 
the scalar field: 2^/^a^/^(y9 = Ph'PL 



3 



is sufficient to lead to the formation of the layered phase through a phase transition. In 
our context we consider this layered phase on the lattice as a possible paradigm on how a 
localisation of the fields, obeying to non-perturbative interactions, may be carried out on 
the brane due to confining interactions in the bulk. 

In this paper we study a simplified realisation of the lattice action given by Eq. (|1.9|) 
below. This is inspired by Eq. fll.7p i.e. to set the fifth (transverse) direction couplings to 
a strong coupling regime while we neglect the explicit role of the warp factor in the lattice 
action. Therefore the lattice action (which leads to the five dimensional Higgs model in the 
naive continuum limit ([Zj, |H]) ) reads in standard notation: 

Sl = 5'gauge + 'S's^alar 

= E (l-cost/^,(a;)) + ^ 5^ /?;(l-cos[/^5(a:)) 

X l</i<j/<4 X l</^<4 

X 

+ Y}{1 - 2(3n - 4/3, - (3'^)^1{x)^l{x) + f3n{^l{x)^L{x))\ (1.9) 

X 

Apart from the resulting simplicity in the context of the phase diagram analysis, a connection 
of this work with previous studies of the layered phase can be achieved. Moreover our 
impression is that the full lattice model is likely to produce physically similar results with 
the present simplified version. This was also the case for the pure U(l) gauge model. The 
'static' representation of the model for which the gauge couplings were fixed by hand gave 
equivalent resuts with the model in which the warp factor was used for the scaling of the 
gauge couplings |4|. 



2 The order parameters and the choice of coupUngs 

We study the abelian Higgs model on the lattice by using numerical methods. The action is 
given explicitly by Eq. (jl.9p . We define five order parameters, making also the distinction 
between space-like and transverse-like ones. These are the following: 

Space-like Plaquette: Ps =< —rT^Y^ cosUu^lx) > 

Transverse-like Plaquette: Pt =< "^^tK cosf/^5(x) > 

X l<At<4 
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Space-like Link: Ls =< ^ ^ cos(x(a; + fi) + Af,{x) - xix)) > 

X l<At<4 

Transverse-like Link: =< cos(x(a; + 5) + ^§(0;) — x{^)) > 

X 

Higgs field measure squared: = ^^p^(x) 

X 

We have assumed the polar form for the scalar field, i.e. (^l = p{x)e^^^^\ 

In jS] this model has been already studied and a first exploration for the phase diagram 
is available. In that work, since the parameter space is very large, consisting of five lattice 
parameters,the choice has been made to fix f3g to 0.5, to 0.001 and consider two values 
of Pr (0.1 and 0.01) and explore the parameter space Ph)- Under these conditions the 
analysis of the order parameters defined above yielded a phase diagram consisting of the 
three expected phases which are the confining phase (S), the Coulomb phase (C5) and the 
Higgs phase (i^s) each of them defined in five dimensions. In addition a fourth phase is 
present: a Higgs phase in four dimensions {H4) (see Fig^. The distinction between H4 and 

can be achieved due to the different behaviour of the transverse-like order parameters 
within the two phases. Details and conclusions on the existence of this layer Higgs phase can 
be found in [Sj. Let us refer also that the identification for the order of the phase transitions 
was possible and has lead to the conclusion that (for the two values of jSn used) both H4, 

are separated from the confining phase by a first order phase transition. We reproduce 
the phase diagram for /^/j = 0.1 as it was depicted in (FigH)). 



3 Searching for a second order phase transition 

At this point the question arises whether it could be possible for the layer phase to 
appear via a second order phase transition. Following 0, we consider the system being in 
the confining regime by setting (3g = 0.5 and fixing to the very small value 0.001. We 
expect the phase transitions to the Higgs phases to be weaker as the Higgs self coupling 
increases. We explore the order of the S — phase transition by setting the transverse 
gauge coupling to 0.2 while we increase Pr. In advance, it should be noted that, as we 
move to larger values of Pr, the relative positions of the phases in the phase diagram are 
substantially similar to what is shown in Fig|T]for Pr = 0.1. So, setting P'g to 0.2, we always 
explore the S — phase transition. 
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Figure 1: Phase diagram of the 5-D Abehan Higgs Model with the space-hke gauge couphng 
set to the strong couphng Pg = 0.5 (taken from (Hj). 



In the sequel we give strong evidence that, the S — first order phase transition line 
ends at a second order point followed by a crossover region. At the same moment the S — 
phase transition remains first order. This additional fact confirms the special nature of the 
four dimensional layer Higgs phase. 

We give now information for the simulating process. We used a 4-hit metropolis algo- 
rithm for the updating of the fields. In addition we implemented the global radial algorithm 
and the overrelaxation algorithm for the updating of the Higgs field. We used four lattice 
volumes, 8^, 10^, 12^, 14^, and we performed 20000-30000 measurements for each point 
which we analyzed in the parameter space. We studied a large number of jSji values before 
concentrating our study to the interval [0.140, 0.165] in which the first order phase transition 
turns to be a weaker one before it passes to the crossover region. 

In the subsequent paragraphs we present our results which are based upon using the 
hysteresis loop technique, the finite volume size scaling, the susceptibility and the study of 
the correlation functions for the Higgs field measure squared. 

3.1 Hysteresis loop technique results 

The first tool for the exploration of the phase diagram with (3fi is the hysteresis loop tech- 
nique. Although this technique gives results that have to be taken into account with caution 
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quantitatively, nevertheless they prove to be very useful as qualitative ones. To this end we 
use the hysteresis loop results as a general guide to get a crude estimate on the Pr interval 
within which the phase transition is converted from a first order to a higher order one. In 
Fig|21we depict the hysteresis loop results for the four -dimensional gauge invariant quantity 
Ls and for four values of Pr, namely Pr = 0.143,0.149,0.153,0.160. The lattice volume 
in this example is 8^. One can see from the figure that while there is a well formed loop 
for Pr = 0.143 indicating a first order phase transition, this changes to a smaller one for 
Pr = 0.149, and it seems to disappear for Pr = 0.153. Although this value should not be 
taken too seriously, one should keep in mind that around the value Pr = 0.153 a weaker 
phase transition is still present. Furthermore we have to mention that the transverse link 
quantity, Ljn, (not shown in the figure) remains almost unaffected by the phase transition, 
being stuck to a very small value close to zero (for details see |8j). 



Pg = 0.20 




0.35 I 1 1 1 1 ' 1 

0.48 0.485 0.49 0.495 0.5 0.505 0.51 



Figure 2: Hysteresis loops showing that the loop for the link-space order parameter disap- 
pears for Pr values bigger than 0.153. 

In Fig. El we give an example of the different phase transition orders of the S — 

and S* — if 5 transitions, both for Pr = 0.158 and lattice volume 8^. In Fig|2t' we present 

the hysteresis loop results on Ps and Pt for P'g = 0.20. The behaviour of Ps indicates 

a phase transition though a smooth one since there is no hysteresis loop, while the Pt is 

almost constant and equals 0.1, in accord with the strong coupling prediction Pg/2. This 

figure should be compared with the Figint", which refers to P'g = 0.80. The hysteresis loop 

results shows a very strong first order phase transition, exhibited by both Ps, Pt^^. This 

M Notice that the unbroken phase is a confining one due to the fact that Ps and Pt follow the strong 
coupling limits f3g/2 and /3g/2 respectively (for more on that see 0). 
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behaviour refers to the S — phase transition. Figures and Eli show the behaviour 
for B? which illustrates the fact that for both cases the system passes to a broken phase. 
In other words, increasing (3 J one finds two different Higgs phases (see for example FigCJ, 
a four-dimensional and a five-dimensional one both separated from the five-dimensional 
confining phase by phase transitions of different orders. 




Figure 3: Characteristic examples showing the obviously different order of the phase tran- 
sition for the cases S — H4 and 5 — if 5 at the same value of Pr, for two different values of 



3.2 Finite volume size scaling 

As it has been discussed in jH] one of the main features of the S — H4 phase transition 
is the multi-layer structure. This means that since the system undergoes a transition to 
a four-dimensional phase rather than a five dimensional one, some special signal should 
appear. Besides a first order phase transition this consists of a multipeak structure in the 
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finite lattice volume histograms for tlie gauge invariant observables, instead of the expected 
behaviour of the two-peak structure. Furthermore it has been shown that every space- 
hke gauge invariant quantity defined on each space-hke volume (i.e. a four-dimensional 
layer) 'feels' the phase transition for different pseudocritical values of the lattice parameters. 
Since this is a consequence of the finite lattice volume used for Monte-Carlo simulations in 
combination with the four-dimensional dynamics when the layer phase arises, we justify the 
choice of analysing the results on the four-dimensional subspace. 

In FiglHwe depict the histograms of the Higgs field measure squared, for = 0.20 
and three values of (3r. All the three histograms refer to (3h values in the critical region. The 
lattice volume in this figure is 14^. The histograms refer to four-dimensional (space-like) 
volume. The two peak structure is more pronounced for the smaller value of (i.e. 0.153), 
where the two peaks are totally separated. For Pji = 0.155 the two peak structure is less 
emphasised while for (3ji = 0.158 it has already disappeared. In order for someone to use this 
method with more safety the lattice volume dependence of the two peak structure should 
be taken into account. This is provided in FigjSl . In Fig|3^ it is easily seen that the two 
peaks become well separated as the lattice length increases from 10 to 14 which serves as 
an indication of a first order phase transition for the case of (3r = 0.153. This has to be 
compared with the really inversed behaviour for = 0.158 shown in Fig|3t. The jS^i = 0.155 
case, Figl^j, for which the peak separation does not change significantly as the lattice length 
goes from 10 to 14, gives an estimate of a first order phase transition becoming much weaker 
and probably of higher order. 
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Figure 4: Histograms of R^ over a space-like volume for three values of (3ii = 
0.153, 0.155, 0.158 and lattice length N = 14. 

Let us now present more quantitative results by giving the results for the susceptibility 
of R^ on the layers for various values of (3^- This is defined by: 
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where Vg denotes the space-hke lattice volume. The results are depicted in FiglHl The 
errors have been calculated by using the Jackknife method. It is known that a first order 
phase transition is signalled by a linear increase of the maximum of the susceptibility with 
the volume. This is actually the case for (3r = 0.149 and 0.153. The situation changes for 
= 0.155 where the linear behaviour is apparently absent. In addition, for the bigger values 
= 0.158 and 0.160 there is not a clear increase with the volume. This case corresponds 
to a crossover behaviour. Therefore, the conclusion is that in the vicinity of (3r = 0.155 
we meet with the well known situation, where a first order phase transition line ends to a 
second order phase transition point followed by a crossover. 
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Figure 6: The susceptibility versus the space-hke volume for five values of in the critical 
/?R region. 



3.3 Correlation functions 

In this section we present the behaviour of two correlation functions, one defined on the 
whole five-dimensional space and the other on the space-like, four-dimensional one. These 
correlation functions involve the Higgs field measure squared i?^, defined in section 2. The 
definition of the correlation functions is given by: 

^^■^^^^ - \ < {R-)j >-< {R-\ >2 ^^-^^ 

where n takes values from 1 to (i.e. the lattice size). The indices S and T are used to 
distinguish the correlators. The one defined in the transverse direction is noted with the 
index T . The other defined in the space-like volume is denoted with S. 

The results for the two correlators are radically different. An example of our results is 
shown in FigjTI This refers to the case of = 14 lattice size for three values of 13^. We see 
that while Ct decreases very fast, reaching zero and fluctuating around it, Cs takes values 
different from zero. This serves as a clear evidence that a layered phase is formed. The layers 
are decoupled as a consequence of the strong coupling imposed on the transverse direction, 
which has the implication of vanishing Ct- Moreover the rather reasonable behaviour of C5 
shows that inside the layers a four dimensional dynamics is still met as it might be expected. 

Another very interesting feature of the Cs correlation function is that as the (Sr value 
decreases the curve becomes more flat. We should note that in the case of a second order 
phase transition and for inflnite volume this should be really flat. This is a fact corresponding 
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Figure 7: The space-like and time-like correlation functions for L = 14 and for three values 
of = 0.155, 0.158, 0.160, in the region of (i^ where the susceptibilities show a peak. 



Lattice size = 0.155 I3r = 0.158 I3r = 0.160 
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0.145(3) 


0.165(4) 


0.181(4) 


12 


0.112(3) 


0.137(3) 


0.167(7) 


14 


0.090(3) 


0.131(6) 


0.152(7) 



Table 1: The masses in lattice units. We observe that for L = 14 and (3r = 0.155 the value 
for the mass parameter has decreased by a factor of 1.7 in comparison with the (3r = 0.160 
corresponding value. 

to infinite correlation length or vanishing mass for the lightest scalar mode. In other words, 
by adjusting the f3h value into the critical region we might expect a mass behaviour of the 
type rris oc {Pr — jS'^y . The light scalar mass calculation can be achieved by using a fit of the 
form const x cosh(ms(x — N/2)) to the correlation functions Cs- The parameter rris is the 
dimensionless mass parameter of the scalar mode. An example of the fits is shown in Fig|7| 
The results for rris for the cases considered are shown in Tabled From that Table and for the 
largest lattice size used we can see that decreases by a factor of 1.7 between /Sr = 0.160 
and 0.155. A more clear signal for the vanishing rris would require bigger volumes and still 
higher computer time. Nevertheless, after considering the previous analysis on susceptibility 
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combined with the results from the study of the correlations, we are justified to estimate 
that at [3r = 0.155(2) a second order phase transition point should be expected. 

4 Conclusions 

We believe that we have serious evidence that the five-dimensional Abelian Higgs model with 
strong coupled interactions along the fifth (transverse) direction reveals a four-dimensional 
dynamics with broken gauge symmetry. This occurs via a second order phase transition. 
The existence of the layered phase can be considered as a realisation for the localisation of 
the gauge and scalar fields for models defined in a higher dimensional space with the extra 
dimensions being warped. Although the lattice volumes and the computer power available is 
not conclusive for the second order critical point (so that the calculation of critical exponents 
is out of consideration for the moment), our results provide an estimate for the value of the 
Higgs self coupling at which the line of the first order transition line ends in a second order 
transition point along the four-dimensional space. 
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